Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 9/08/03 



O 
O 



(N 
> 
(N 

m 

(N 

m 
O 

o 

a: 
6 

c3 



13 



THE SPEED OF THE 'BULLET' IN THE MERGING GALAXY CLUSTER 1E0657-56 
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ABSTRACT 

Deep Chandra exposures of the hot galaxy cluster 1E0657-56 have revealed a highly dynamical state of the 
system due to an ongoing merger with a massive subcluster. The system is observed shortly after the first 
core-passage of the infalling subcluster, which moves approximately in the plane of the sky and is preceded 
by a prominent bow shock with Mach number A4 ~ 3. The inferred shock velocity of ~ 4700 kms" 1 has been 
commonly interpreted as the velocity of the 'bullet' subcluster itself. This velocity is unexpectedly high in 
the ACDM cosmology, which may require non-trivial modifications in the dark sector such as additional long- 
range scalar forces if taken at face value. Here we present explicit hydrodynamical toy models of galaxy cluster 
mergers which very well reproduce the observed dynamical state of 1E0657-56 and the mass models inferred 
from gravitational lensing observations. However, despite a shock speed of 4500 kms" 1 , the subcluster's mass 
centroid is moving only with ~ 2600 kms" 1 in the rest frame of the system. The difference arises in part due 
to a gravitationally induced inflow velocity of the gas ahead of the shock towards the bullet, which amounts 
to ~ 11 00 kms" 1 for our assumed 10 : 1 mass ratio of the merger. A second effect is that the shock front 
moves faster than the subcluster itself, enlarging the distance between the subcluster and the bow shock with 
time. We also discuss the expected location of the lensing mass peak relative to the hydrodynamical features 
of the flow, and show that their spatial separation depends sensitively on the relative concentrations and gas 
fractions of the merging clusters, in addition to being highly time dependent. A generic ACDM collision 
model, where a bullet subcluster with concentration c = 7.2 merges with a parent cluster with concentration 
c = 3 on a zero-energy orbit, reproduces all the main observational features seen in 1E0657-56 with good 
accuracy, suggesting that 1E0657-56 is well in line with expectations from standard cosmological models. 
In theories with an additional "fifth" force in the dark sector, the bullet subcluster can be accelerated beyond 
the velocity reached in ACDM, and the spatial offset between the X-ray peak and the mass centroid of the 
subcluster can be significantly enlarged. Our results stress the need for explicit hydrodynamical models for the 
proper interpretation of actively merging systems such as 1E0657-56. 

Subject headings: clusters: structure - clusters: x-ray observations - methods: numerical 



1. INTRODUCTION 

The massive cluster of galaxies 1E0657-56 at z = 0.296 
offers a so-far unique setting for testing many aspects of 
our understanding of dark matter and baryonic physics. In 
this system, a massive sub-cluster (the "bullet" with M^oo — 
1.5 x 1O 14 M ) has fallen through the main cluster (M200 — 
1.5 x 10 15 M^) on a trajectory n e arly exactly in the p lane of 
the sky (IMarkevitch et alJ l2002t iBarrena et afll2002l) . The 
matter distribution and masses quoted above h ave been de- 
termined from weak and stron g lensing studies dClowe et alJ 
12004 12006t iBradac et alj|2000) . Fitting the large-field weak 
lensing data with two spherically symmetric NFW mass dis- 
tributions leads to (D. Clowe, private communication) r2oo = 
2136 kpc, c = 1.94 and r 200 = 995 kpc, c = 7.12, for the 
main and subcluster respectively, taking ho = 0.7, f2 m = 0.3, 
fl\ = 0.3. The distribution of baryons in the diffuse intra- 
cluster medium ( ICM) is obtained from 500 ks Chandra X- 
ray observations (Markevitch 2006). By fitting the disconti- 
nuity in the gas properties across the prominent bow shock 
preceding the gas bullet, the shock Mach number has been 
determined to be M = 3.0 ± 0.4 (IMarkevitch et alJ 120021: 
Markevitch 2006), which together with a measurement of the 
pre-shock temperature tra nslates to an infe rred shock velocity 
of v s = 4740!li!!kms" 1 (Markevitch 2006). Note that the er- 



rors quoted in (Markevitch 200jl) were symmetrized and the 
central velocity rounded down; the values used here are the 
most accurate ones (M. Markevitch, private communication). 

In previous work, it has been assumed that this shock ve- 
locity is equal to the su bcluster's relative velocity with re- 
spect t o the parent cluster ^Markevitch et al.l2 002: Mar kevitch] 
I2006t lHayashi &Whitd 120061 among others'). This veloc- 
ity is much higher t han naturally expected according to the 
ACD M cosmology dHavashi & White! 120061: iFarrar & Rosenl 
|2006|). Using the larg e cosmological Mill e nnium simulation 
(ISpringel et al.1 l2005bl) . lHayashi & White! (120061) have mea- 
sured the distribution function of the velocity of the most mas- 
sive dark matter substructure in massive clusters of galaxies, 
fitting the results with 
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where V200 is the circular velocity of the parent at the 
virial radius, and /(> V su b) gives the fraction of clus- 
ters with a most massive substructure moving faster than 
V su b- Adopting for V su b a shock velocity of 4500 kms" 1 
(IMarkevitch et al.ll2004l) and for the virial velocity of the par- 
ent cluster V200 = 2380km s" 1 (based on a mass model with 
M 20 o = 2.16x lO' S /r'Af p), R 2 ao = 1.64/r'Mpc and h = 0.73), 
lHavashi & White! d2006l) found a relative likelihood of / ~ 
0.01 for 1E0657-56, suggesting that the system, while rare, 

can be plausibly obtained in ACDM. 

How ever, as previously pointed out by Farrar & Rosen 
(2006) this estimate appears overly optimistic in light of 
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the new improved data on 1E0657-56. This data has in- 
creased the shock velocity to 4740 km s" 1 , while the mass- 
estimate for M200 has seen a downward revision by a fac- 
tor of ~ 2. The virial velo city V200 that is no rmally as- 
sociated with a mass M200 (Navarro et al.l fl996h is given 
by V200 = [10Gi/(z)M 2 oo] I//3 - Accounting for the revised 
lower mass estimates from gravitational lensing dClowe et al.l 
2006; iBradac et al.ll2006l) . V200 goes down to at most V200 ~ 
1770kms _1 at the cluster redshift of z = 0.296. This means 
that the critical ratio V su b/V200 changes from 1.9 to 2.68, im- 
plying that the estimated likelihood of seeing a subcluster as 
fast as the bullet in 1E0657-56 is really only / ~ 0.8 x 10" 6 , 
based on Equation (HJ. Given the low abundance of clusters as 
massive as 1E0657-56, this is an exceedingly small probabil- 
ity which cannot be straightforwardly understood in ACDM. 

An attempt to estimate the relative velocity of the sub- 
cluster and main cluster when separated by 720 kpc, taking 
reasonable initial conditions and the measured mass distribu- 
tions, finds a maximu m relative velocity of about 3500km s" 1 
dFarrar & Rosenll2006l) . suggesting the possible presence of 
an additional long-range force between dark matter con- 
centrations if 4740 k ms" 1 is the true relative velocity (see 
Farrar & Rosen 2006, and references therein). The results of 
Havashi & White (2006) and the conclusion to be drawn from 
Farrar & R osen (2006) depend critically on the identification 
of the velocity of the dark-matter bullet with the observation- 
ally inferred shock velocity. Given the importance of a reli- 
able determination of the bullet's velocity, we here explore the 
relation between the true velocity of the dark matter subcluster 
and the properties of the shock in the gas distribution, using 
simplified 'toy' models of the collision to make the physics 
clear. 

We briefly summarize our methodology in $2] An anal- 
yses of the hydrodynamical state of a best-matching merger 
simulation is given in §|3] which illustrates in particular that 
a significant difference between the shock and bullet veloc- 
ity is naturally expected. In §|4] we consider slight variations 
of our default model and quantify in more detail the quality 
of the match to the observed 1E0657-56 system. This is fol- 
lowed with a discussion of the dependence on further model 
parameters in ij5] In ^6]we examine simulations that include a 
fifth force between dark matter particles. Finally, we conclude 
with a discussion in $7] 

2. NUMERICAL TOY MODELS FOR THE BULLET CLUSTER 

As our default model, we consider the collision of two iso- 
lated galaxy clusters at z = of mass M 20 o= 1.5 x 1O 15 M 
andM2oo = 1.5 x 10 14 M Q , respectively, both containing a uni- 
versal baryon fraction of fi, = ttb/tt m = 0.17, consistent with 
the lat est cosmological co nstraints from the WMAP 3rd-year 
data dSpergel et al.l 120061) . We take the clusters to be spher- 
ically s ymmetric, with th e dark matter following an NFW- 
profile dNavarro et al.l 19971) . truncated at the virial radius, and 
the gas being in hydrostatic equilibrium with a small core 
in the center, but otherwise following the NFW profile as 
well. For the concentrations, we adopt the values c = 2.0 and 
c = 7.2 suggested by the mass fits of Clowe (2006, private 
communication). Assuming a conventional virial overden- 
sity of 200 with respect to the critical density, and a Hubble 
constant of Ho = 73kms~ 1 Mpc~ 1 , the virial velocity and ra- 
dius of the parent cluster are 1680 kms" 1 and 2300 kpc, while 
those of the subcluster are 780kms _1 and 1070kpc, respec- 
tively. Numerical realizations of these equilibrium clusters as 
N-body/SPH models are constructed with the techniques de- 



scribed in Spri ngel et al.l d2005al) . 

We assume that the clusters fall together from infinity on 
a zero-energy orbit, giving the infalling subcluster a veloc- 
ity of 1870 kms -1 in the center-of-mass frame, at the time 
when the virial radii first touch. This is also the time when 
we start our simulations. Note that all the absolute velocities 
we quote are relative to the center-of-mass frame, hence the 
parent cluster moves with -187kms _1 at the start. The sim- 
ulations themselves are evolved with a novel version of the 
cosmologic al smoothed particle hydrodynamic s (SPH) code 
GADGET2 dSpringel et al.ll200U ISpringel 2005b- In our stan- 
dard simulations, we employ 4.4 million particles, equally di- 
vided between dark matter and gas, and with equal mass reso- 
lution for the incoming subcluster and the parent cluster. This 
results in a gas mass resolution of m gas = 1 .275 x 10 8 M Q , and 
a dark matter resolution of m^m = 6.225 x 10 8 M Q . We adopt 
a gravitational softening length of 2.8 kpc, which roughly de- 
marks our spatial resolution limit. 

In our default calculation we consider a direct head-on en- 
counter of the two clusters, making the problem effectively 
axially symmetric around the collision axis. The need to ac- 
count for the non-linear response of the dark matter during 
the merger requires three-dimensional simulations neverthe- 
less. We have also considered a number of simulations with 
a small amount of orbital angular momentum in the collision, 
measured by the minimum separation b the clusters would at- 
tain at perihelion if they were point masses. Also, we con- 
sidered variations of our initial conditions where the relative 
concentrations and gas fractions of the initial cluster models 
were altered, or where a fifth force between dark matter parti- 
cles was included. Finally, for testing numerical convergence, 
we have also carried out selected simulations with 35.2 mil- 
lion particles, resulting in 8 times better mass resolution, and 
a 2 times better spatial resolution in each dimension. 

Since the measured temperature profile across the shock 
is consistent w ith instant shock-heating of the electrons 
(Markevitch 2006), we model the ICM as an ideal gas with an 
adiabatic index of 7 = 5/3. For simplicity, we neglect radia- 
tive cooling of the gas, and consequently also star formation 
and the stellar systems of the cluster galaxies. We also disre- 
gard magnetic fields, which could in principle play an impor- 
tant role in stabilizing cold front against fluid instabilities and 
thermal conduction (lEttori & Fabian] l2000t IVikhlinin et al.l 
1200 lbt iLvu tikov 2006). Our simulations are carried out in 
physical coordinates, ignoring the cosmological background 
expansion, the infall of other structures, and the accretion of 
diffuse matter. While our approach is much more restricted 
than fully self-consistent cosmological modeling, the short 
time required for completing the merger (of the order of a 
few Gyrs) should make it a reasonably good approximation. 
At the same time, the simplicity of the toy model makes the 
physical phenomena easier to identify, in particular with re- 
spect to the relation between observed shock properties and 
the underlying kinematics of the dark matter, which is our 
main interest here. We shall begin with the simplest of the toy 
models, and then explore how the conclusions change as the 
model is varied. 



3. RESULTS FOR THE BASIC TOY MODEL 

In Figure Q] we show the time evolution of our default sim- 
ulation model in the form of maps of X-ray surface brightness 




FIG. 1. — Time evolution of our default merger simulation, as seen in X-ray surface brightness (top) and in the luminosity-weighted projected temperature 
(bottom). The panel in the middle of both series of panels shows a time close to the best-match to the observed cluster 1E0657-56. 
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FIG. 2. — X-ray surface brightness of 1E0657-56 observed with Chandra (left). Clearly visible are two X-ray peaks. The wedge-like structure on the right 
associated with the 'bullet' is bounded by a sharp contact discontinuity, which is interpreted as a cold front. The prominent bow shock in front demonstrates that 
the subcluster is moving to the right with high velocity. The panel on the right-hand side shows the X-ray surface brightness in one of our merger simulations, 
roughly drawn on the same scale and with similar dynamic range in the color table. 



and projected luminosity- weighted temperature 3 . We only 
consider simple bremsstrahlung, with an X-ray volume emis- 
sivity ex oc PgVT, where p g and T are the local gas density 
and temperature, respectively. The bullet cluster enters the 
parent at time t = from the left, and passes the core at around 
t 1200Myr. The bullet is preceeded by a prominent bow 
shock, which develops a 'bend' during core passage, when 
the central region of the shock front speeds up. The compres- 
sion of the central gas during core passage makes the centers 
of the parent and the bullet cluster very bright in X-rays for a 
few hundred Myr s. In addition, the bull et develops a promi- 
nent 'cold front' (Vikhli nTn et al.ll2001al) after core passage, 
which can be seen as a sharp wedge-like 'edge' in the X-ray 
surface brightness. The temperature maps show the strong 
heating of the gas by the shock front, while the bullet itself 
is seen to be quite cold, revealing th at the wedge-like edge of 
the bu llet is really a cold front (see iMarkevitch & Vik hlinin 
120071 for a review). Note that the maps shown in Fig. [Hex- 
tend to substantially lower surface brightness levels than are 
observationally accessible. 

In Figure[2] we show an observed X-ray ma p of 1E0657-56 
(based on data presented in Markevitch 2006) and compare it 
to a map of the simulation drawn with a dynamic range that 
roughly matches that accessible in the observations. We se- 
lected a time that provides approximately the best match to 
the observed X-ray morphology, as will be discussed in more 
detail later. The simulation clearly shows a double peaked 
structure in the X-ray emission similar to that observed, with 
a wedge-like shape for the bullet that moves to the right, and 
with a prominent bow shock in front of the bullet. This broad 
agreement reconfirms the interpretation of 1E0657-56 as a 
merging system. We now turn to a more detailed analysis of 
the simulation results, in order to evaluate how well they can 
quantitatively reproduce the properties of the observed sys- 
tem. 



3.1. The dynamical structure of the merging system 

3 A video of the merger simulation is available for download at 

http : / / www . mpa-garching . mpg . de/~volker /bullet 



A much clearer understanding of the hydrodynamics in- 
volved is obtained by looking at profiles of the gas and the 
dark matter distribution along a line through the nose of the 
shock. This is done in Figure [3] for the local X-ray emis- 
sivity, gas velocity field, temperature, pressure, gas density, 
and total projected mass. In this figure at time t = 1340Myr 
(corresponding to the central panel in the time evolution of 
Fig. 1), we can clearly identify the location of the bow shock 
at x s ~ 580kpc, and a contact discontinuity at x c ~ 482kpc. 
The latter marks the edge of the bullet; behind it, the gas is 
quite dense and cold, and very bright in X-rays. The pressure 
is essentially continuous across this front. 

At the bow shock, the density, velocity, pressure, temper- 
ature, and X-ray emissivity are all discontinuous (as is the 
specific entropy). The sizes of the corresponding jumps are 
consistently described by the Rankine-Hugoniot relations for 
an adopted Mach number of M. = 2.8, as indicated in the pan- 
els of the figure. The preshock-temperature is ~ 9kev, just 
what is observed, corresponding to a pre-shock soundspeed 
of c s = 1600kms _I . The inferred shock velocity is hence 
v a = Mc s = 4488 km s" 1 . Note that a slightly higher Mach 
number of 3.0 matches the numerical results still reasonably 
well, and would then in fact result in a shock velocity exactly 
equal to the value inferred from the most recent observational 
estimates (references). In any case, the value of M. = 2.8 is 
well inside the observational uncertainty. 

We note that the opening angle of the Mach cone seen in 
our simulations is much wider than expected in the simple 
situation of a perturber moving with constant supersonic ve- 
locity in a homogeneous medium that is at rest. This is in part 
a result of the velocity structure of the upstream gas, which is 
not at rest, a point we shall return to later on. In addition, the 
speed of the bullet is not constant, and the temperature field is 
not strictly homogeneous. As a combined result of these com- 
plications, the opening angle of the bow shock cannot easily 
be used to obtain an independent measure of the strength of 
the shock. 

3.2. The velocity of the bullet 

An important point to note about Figure [3] is that the shock 
velocity itself does not give the expected jump in the velocity 
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FIG. 3. — Local profiles of different hydrodynamical quantities along the symmetry axis through the nose of the shock, at time t = 1340Myr after start of 
the simulation, corresponding to 1 80 Myr after core passage. From top to bottom, we show the X-ray emissivity, gas streaming velocity along the ^-direction, 
gas temperature, gas pressure, and gas density. Finally, the bottom right panel shows the projected total mass surface density along the ,t-direction. The dashed 
vertical line indicates the location of the bow shock, while the horizontal dashed lines give the expected jump in the respective quantity for a Mach number of 
M = 2.8. The vertical dotted line marks the location of the edge of the bullet. 
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FIG. 4. — The top panel shows the location of the bow shock (red) and the 
edge of the bullet (blue) as a function of time. Also shown are the positions 
of the mass centroids of the bullet and the cluster with solid lines (as labeled). 
In the bottom, the corresponding velocities are shown, obtained by numerical 
differentiation of the measured locations as a function of time. The time of 
core passage is marked with a dashed vertical line, while the time of the best 
fit of our simulation to the observed bullet cluster is marked with a dotted 
line. 



of the gas across the shock. The latter is only 



Av = Vi - V2 = v s 



1 Pi I 

1 I = v.. 

Pi 



2(7W 2 -1) 
( 7 +l)A^ 2 



= 2937 km s" 



(2) 

as a result of the continuity equation. Here p\ and p2 give 
the gas densities in front of and behind the shock, and vi and 
V2 are the corresponding velocities. If the latter are measured 
in the shock's rest frame, we have p\ v\ = P2V2, which leads 
together with the Rankine-Hugoniot jump condition for p\/ p2 
to Equation (f2j). 

But perhaps the most noteworthy feature in Figure [5] is the 
fact that the velocity of the gaseous bullet, which is located 
at around ~ 400 kpc, is not at all equal to the shock veloc- 
ity v s , as has often been assumed. Furthermore, our results 
show that the mass centroid of the subcluster is in fact mov- 
ing much slower; its velocity lies around ~ 2600km s" 1 , yet 



the inferred shock speed is v s = 4488 km s The difference 
originates primarily for two reasons: 

1 . The pre-shock gas is not at rest. Instead, it falls towards 
the bullet, reaching nearly Vi n f a u = llOOkms" 1 when it 
encounters the shock. This is because the gravity of 
the incoming bullet is felt already ahead of the shock 
front, unlike hydrodynamical forces. We shall discuss 
this point further below. For now we note that this 
means that the shock front moves only with a velocity 
Vj = Vj - Vinfaii in the rest-frame of the system. 

2. The subcluster's mass centroid and the X-ray centroid 
are not moving as fast as the shock front itself, and are, 
in fact, increasingly falling behind. Note that the post- 
shock gas has only a velocity of Av relative to the pre- 
shock gas, and not the full speed v s of the shock. The 
bullet will hence experience significant braking by ram- 
pressure if it moves faster than Av - Vin&n in the rest- 
frame, which it does. This suggests that a difference 
between the speed of the bullet and the shock front may 
naturally arise, but the size of this effect is hard to es- 
timate due to the time variability and the complicated 
geometry of the merger. But as we quantify below, a 
velocity difference between the subcluster's mass cen- 
troid and the bow shock is indeed present and amounts 
to ~ 700 km s" 1 at the time of best match. 

In Figure |4j we show the locations of the shock front, the 
edge of the bullet, and the mass centroids of parent and sub- 
cluster as a function of time, from a brief moment before 
core passage over a period of 400 Myr. Clearly seen is that 
the shock front moves faster than the mass centroid asso- 
ciated with the bullet. In fact, numerical differentiation of 
these curves shows that the shock front is moving with v s ~ 
3300 kms -1 in the center-of-mass frame at time t = 1340 Myr, 
while the mass centroid is moving with 2600km s" 1 . The edge 
of the bullet stays closer to the shock front, however, only ini- 
tially their distance grows to ~ 80 kpc but then stays roughly 
constant. Given the pre-shock infall of 1 100 kms -1 of the gas, 
the velocity of the shock front in the rest-frame of the up- 
stream gas is ~ 4400 kms" 1 , consistent with the shock veloc- 
ity v s = Aic s inferred from the Rankine-Hugoniot jump con- 
ditions. 

We can also see that the core passage has happened 180 Myr 
before the time displayed in Figure[5J at which point the mass 
peaks have a separation of 627 kpc. This is a bit smaller than 
the 720 kpc cited for the separation of the mass centroids of 
1E0657-56, based on gravitational lensing analysis. How- 
ever, just 60 Myr later, the distance grows to this level, with 
an unchanged shock strength, so that the time of best match 
using this measure could be defined as t ~ 1400 Myr. 

Comparing with the positions of the hydrodynamical fea- 
tures in Figure [4] we see that for the chosen simulation pa- 
rameters the mass peak of the bullet lies slightly ahead of 
its X-ray edge at the time of best match, at which point it 
is in fact also quite close to the location of the shock front. 
While this is broadly consistent with the results from weak- 
and s t rong lensing analys is of 1E0657-56 (IClowe et alJl2004l 
2006; Bradac et al. 20061), the size of the offset in this simula- 
tion is clearly smaller than inferred observationally. We note 
that this is also the case for the offset between the X-ray emis- 
sion and mass peak associated with the parent cluster: those 
features are quite well aligned in our simulation model, un- 
like in the observations. However, we will consider a small 
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FIG. 5. — The gravitational potential along the symmetry axis of the merger, shown at different times. In each panel, the red line (top-most line) shows the 
potential generated by the mass belonging initially to the infalling sub-cluster, while the blue line gives the potential generated by the parent clusters mass. The 
black line shows their sum (bottom-most line). The panel at t = 1340 Myr additionally contains the potentials generated by the original mass distributions of 
bullet and parent (dashed lines), horizontally translated to their current positions. The difference in the two total potentials arises from rearrangements in the 
mass distributions of the two components. 



variation of our default simulation model later on in $4]where 
a more accurate match to the observational data is obtained. 

3.3. Gas infall towards the bullet 

In the previous section we have seen that the upstream gas 
ahead of the bullet's shock front is not at rest relative to the 
main parent cluster. Rather it is infalling towards the in- 
coming bullet with a significant velocity. This effect con- 
tributes substantially to the large difference between the in- 
ferred shock velocity and the actual velocity of the infalling 
subcluster. Here want to examine more closely how this ve- 
locity arises. 

In Figure |5j we show the time evolution of the gravita- 
tional potential along the symmetry axes during the merger. 
Each panel shows the potential at a different time, both for 
the total mass and separately for the material that originally 
defines each of the two clusters before they start overlap- 
ping. In the early phases of the collision, the mass distri- 
bution of each of the clusters changes relatively little, such 
that the total potential is approximately given by translating 
the initial potentials of the two components, and summing 
them up. However, during core passage the summed poten- 
tial changes rapidly in time, leading to a phase of violent re- 
laxation that reshuffles the energies of individual dark matter 
particles. This irreversibly modifies the mass distributions of 
the involved components, and hence their potentials. Recall 



that in a static gravitational potential the specific energies of 
individual dark matter particles, e, = v?/2 +<£>(%•), are constant 
along their orbits Xi(t), but if the potential varies in time, their 
energies evolve at a rate de,/df = d^/dt, i.e. the particle ener- 
gies simply change with the rate of change of the potential at 
their current location. It is precisely during core passage when 
the potential fluctuates on a very short timescale such that a 
widespread rearrangement of particle energies can occur. 

We now relate this phenomenon to the large infall velocity 
we have found in the gas ahead of the shock. We argue that 
this occurs largely as a response to the change in the gravi- 
tational potential of the system induced during core passage. 
Consider to this end the panel at time 1 340 Myr in Figure[5] In 
addition to the direct measurement of the total potential (solid 
lines), we have also included an estimate of the potential ob- 
tained by linearly translating the potentials of the two original 
mass components to the current positions of the mass cen- 
troids (dashed lines), and then adding them. In the resulting 
comparison, it is seen that the dynamically evolved potential 
has clearly become deeper during core passage relative to that 
expected for unaltered mass distributions in the two compo- 
nents. Note in particular that the total potential in the infall 
region of the gas ahead of the shock front has become deeper 
as well; this implies that the gas in this region is accelerated 
towards the shock considerably more strongly than predicted 
in a model where the dynamical response of the parent and 
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FIG. 6. — Gas infall velocity and mean streaming velocity of the parent 
cluster's dark matter along the symmetry axis of the system. The location of 
the bow shock is clearly seen as vertical discontinuity in the gas velocity. The 
solid black line shows an estimate for the mean dark matter velocity obtained 
from the change in the gravitational potential generated by the time variation 
of the mass distributions of parent cluster and infalling bullet. 



bullet during core passage is neglected. 

We can use the above measurement of the 'excess change' 
A$ in the gravitational potential to obtain a rough estimate 
of the induced dark matter streaming velocity in the parent's 
cluster material in the infall region ahead of the shock. To this 
end we invoke energy conservation for the radial streaming 
velocity of the dark matter, which at a given location should 



roughly increase as v r 



-2A$. This estimate neglects the 



fact that the particles really have a velocity distribution func- 
tion, but should be good for an order of magnitude estimate. 
In Figure [6] we show the result, and compare it with the real 
dark matter streaming velocity measured in the simulation, 
and the mean gas velocity. We see that the above estimate 
matches the measurement quite well. The gas inflow velocity 
is somewhat lower than the induced dark matter streaming. 
This can be understood as a result of gas pressure effects that 
resist the adiabatic compression associated with the induced 
infall, and which are absent for the collisionless dark matter. 

4. IMPROVING THE MATCH TO 1E0657-56 

We now discuss in more detail the quality of the match of 
our default merger model to the observations of 1E0657-56, 
and in particular consider the sensitivity of our results to the 
assumed concentration of the parent cluster. 

In Figure [7] we compare maps of the X-ray emission with 
overlaid contours of the projected total mass for different 
choices for the initial concentration of the parent cluster prior 
to the merger. The value of c = 2.0 adopted in our default 
model (as suggested by new mass models of Clowe matched 
to the lensing data) is quite low for a cluster of this mass, 
given that the typical concentration expected in AC DM mod- 
els for a cluster of this m ass lies around c ^ 4-5 dEke et al.l 
I200U bullock et al. 2001). We therefore calculated a number 
of simulations with higher concentrations, using c = 3,4, and 



5. A larger concentration of the parent system will make it 
more difficult for the gas bullet to penetrate through the center. 
In particular, we expect that the incoming bullet experiences 
a higher ram pressure and will presumably also be stripped 
more efficiently of its dark matter during core passage. 

The results shown in Figure [7] confirm this expectation: 
With increasing concentration, the separation between the X- 
ray bright region of the bullet (where the gas density peaks) 
and its associated dark matter mass peak becomes larger. In 
addition, the shape of the contact discontinuity (the edge of 
the bullet) changes noticeably as well; it becomes progres- 
sively rounder with increasing concentration, such that the 
opening angle of the narrow wedge obtained in the c = 2.0 
case is first enlarged, and finally becomes transformed into a 
round, blunt shape. On the other hand, the shape of the bow 
shock remains essentially unaffected by the change in concen- 
tration. 

By visual inspection, it appears that the c = 3.0 case among 
this set of simulations can best match the observed shape of 
the bullet's edge in 1E0657-56. For this simulation we also 
find a significant spatial offset between the mass peak and 
the X-ray emission, quite simil ar in size to t hat fo und in the 
gravitational lensing analysis of lClowe et al.l(l2006l) . Figure [8] 
quantitatively confirms this. In it we show the location of the 
edge of the bullet, the bow shock, and the two mass peaks, 
as a function of time, both for the default case of c = 2.0 
and the c = 3.0 model. In order to identify the time of best- 
match to the observed system, we plot the positions relative 
to the coordinate of the edge of the bullet and we mark the de- 
sired spatial offsets of the observed system with shaded bands. 
We focus on reproducing the 720 kpc distance on the sky be- 
tw een the mass peaks of parent cluster and subcluster given 
by IClowee t al. (2006), and on the separation of ~ 100 kpc 
between the bow shock and the edge of bullet. For a good 
match, we hence require that these two relative differences 
are reproduced simultaneously. We see that this is quite well 
the case for the c = 3 simulation, at around t = 1420Myr. Our 
default c = 2 simulation underpredicts the separation between 
the bullet's mass peak and edge of the bullet at all times, while 
the c = 4 and c = 5 cases (not shown) overpredict it during the 
times when the mass peaks have the right separation. 

It is interesting that the c = 3 case also reproduces other ex- 
isting observational constraints on 1E0657-56 quite well. In 
Figure|9] we show an enlargement of the simulation's temper- 
ature structure across the shock front and the contact disconti- 
nuity, and compare it to measur ements of the tempe rature pro- 
file of 1E0657-56 obtained by iMarkevitchl (2006) based on 
deep Chandra data. The shock seen in the c = 3.0 simulation 
has a Mach number of M. = 2.9, consistent with the observa- 
tional data. The fall off of the temperature at the cold front 
is reproduced quite well. Also, the pre-shock temperature of 
9 keV is matched well by the model; only further ahead of 
the shock the simulated temperature appears marginally lower 
than the measurements. 

We can also consider the gas masses in spherical apertures 
of 100 kpc around the mass centroids of bullet and parent 
cluster, and about the center of X-ray emission of the bul - 
let. These masses have been measured bv lClowe et alj d2006l) . 
who estimate (2.7 ±0.3) x 1O 12 M of gas around the bullet's 
mass centroid (actually they used the BCG position, which 
we take here as a proxy for the mass centroid), while the 
simulation gives 1.5 x 10 12 Mq. For the parent mass cen- 
troid they get (5.5 ± 0.6) x 1O 12 M , and the simulation yields 
6.1 x 10 12 M Q . Finally, for the bullet's X-ray centroid they 



Fig. 7.— 



X-ray surface brightness and projected mass maps (contours) for simulations with different concentrations of the main cluster, as labeled. 



find (5.8 ± 0.6) x 10 12 M Q while the simulation model gives 
also 5.8 x 10 12 M Q . We note howe ver that unlike foun d in 
the lensing mass re constructions of IClowe et~ai1 d2006l) and 
IBradac et al.l (|2()06), our simulation does not show a signif- 
icant spatial offset between the mass centroid of the parent 
cluster and its associated X-ray centroid. 

In Figure [10] we consider the projected cumulative mass 
profile as a function of radius around the bullet's and the par- 
en t's mass centroids, a nd we compare this to the measurement 
of IBradac et al.l (2006) based on their weak- and strong lens- 
ing reconstruction. Our model matches the result of the lens- 
ing analysis for the parent cluster very well, but predicts sig- 
nificantly less enclosed mass for the bullet itself. This could 
perhaps be remedied by simulations with a more massive, a 
more concentrated, or a non-spherical bullet. Alternatively, it 
could be due to a systematic overestimate in the lensing mass 
estimate of the bullet. While the latter seems unlikely due 
to the extensive tests that have been applied to the mass re- 
construction methods applied to 1E0657-56, the fact that the 
baryonic masses of the simulation model agree well with the 
directly measured gas masses from the X-ray observations is 
a point supporting the simulation's mass model. 

Finally, an interesting observation about the shape of the 
bow shock in our simulated mergers is that it features a clear 
bend towards a larger opening angle, approximately 600 kpc 
away from the symmetry axis of the merger. This feature orig- 
inates during the core passage of the bullet, when the cen- 
tral part of the bullet is accelerated to its peak velocity and 



shoots out on the other side of the core. It is presently unclear 
whether such a feature exists in the observations, as it should 
be located in a region where the available X-ray maps have 
extremely low surface brightness and are dominated by noise. 
With some good will , one can spot a te ntative hint for it in 
the observed maps of Ma rkevitchl $2006), but this is certainly 
far from being significant. Future X-ray observations should 
be able to probe to fainter surface brightness levels, in which 
case this generic prediction of our simulation models could be 
tested. 

4.1. The future evolution of the bullet cluster 

It is interesting to use our simulation model to study the 
expected future evolution of the bullet cluster system. Adopt- 
ing the c = 3 simulation as our best-matching model, we pre- 
dict that the bullet will continue to enlarge its distance from 
the parent cluster for another 900 Myr from the time of best- 
match. After turn around, it falls back and experiences a sec- 
ond core passage 2500 Myr after the first one. This is followed 
by a third passage another 1700 Myr later, after which the bul- 
let is nearly completely destroyed, such that a 'ring down' 
with a quick succession of possibly several further core pas- 
sages leads to a complete coalescence of the cores. This al- 
lows the cluster to reestablish hydrostatic equilibrium, so that 
we can expect 1E0657-56 to reach a fully relaxed state again 
within the next ~ 4Gyr. 

In Figure QT| we show the time evolution of the expected 
bolometric X-ray emission of the system until the merger is 
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FIG. 8. — Relative distances between the mass-centroids of bullet and parent cluster, and the position of the bow shock, relative to the edge of the bullet (i.e. 
the cold front), as a function of time. Adopting a separation of 720 kpc for the distance of the mass centroids, and 100 kpc for the offset of the bullet's mass 
peak relative to the edge of the bullet, we can identify the time of the best match of the simulation to the observations. A simultaneous match of the two target 
separations is obtained for the c = 3 model at time t ~ 1420 Myr (see shaded target regions). The default model with a lower concentration of the parent cluster 
of c = 2.0 provides a much poorer fit, primarily because no significant offset between the position of the bullet's mass and its X-ray edge develops. 



complete. For simplicity, we have simply computed the quan- 
tity L X (x J PgVT dx as a proxy for the bremsstrahlung emis- 
sivity from the gas of the simulations, without taking into ac- 
count metal lines. The vertical dotted lines mark core pas- 
sages of the bullet cluster. These happen slightly earlier than 
the induced peaks in the X-ray emissivity. At the time of best 
match with the observed state of 1E0657-56 (marked with 
a dashed line), the system is in an overheated state, with a 
mean emission-weighted temperature of around 1 1 keV. The 
cluster responds to the compression from the merger by adia- 
batically overexpanding past its equilibrium in the subsequent 
evolution, leading to a comparatively cool state for a period of 
order 1 Gyr before the second collision. Then a further heat- 
ing event occurs when the surviving part of the bullet passes 
through the core a second time. Finally, the third core pas- 
sage leads to near complete disruption of the remains of the 
bullet. We are then left with a relaxed cluster with a mean 
emission-weighted temperature of ~ 7.5 keV. 

5 . DEPENDENCE ON MODEL PARAMETERS 

5.1. Orbital angular momentum of the encounter 

Based on the visual appearance of 1E0657-56 and 
its small line-of-sig ht velocity difference of ~ 600km s" 1 
dBarrena et al .112002 ). it appears that the merger is very nearly 
central, and closely aligned with the plane of the sky. Indeed, 
a slight offset from a central collision can easily introduce a 
significant degree of axial asymmetry larger than observed in 
the cones of bow shock and cold front. We illustrate this in 
Figure [12] where we show in the top panel a projected X-ray 



map for a merger with small orbital angular momentum, pa- 
rameterized by an 'impact parameter' b = 12.5 kpc. Here b 
is the distance at pericenter the clusters would reach if they 
interacted as point masses during the collision. Comparing 
this result with the observed X-ray map in Figure [2] we see 
that the shock structure is already too asymmetric to provide 
a good match to the observed X-ray morphology. 

However, we can always recover symmetry by turning the 
orbital plane from observing it face-on to observing it edge- 
on. In the middle and bottom panels of Figure Q~2]we show the 
same simulation when viewed 45 degrees and 90 degrees with 
respect to the orbital plane. While in these cases symmetry 
is partially/fully recovered, the edge of the bullet and the bow 
shock appear rounder and less sharp, and the edge of the bullet 
develops a shape that agrees comparatively poorly with the 
observed one. The line-of-sight velocity differences obtained 
in the two cases amount to 235 km s" 1 (for 45 degrees) and 
340km s" 1 (for 90 degrees), respectively, which means that 
one will quickly run into tension with the observations when 
one tries to 'hide' a larger orbital angular momentum in the 
encounter than considered here with an edge-on view of the 
system. 

We hence conclude that 1E0657-56 should be quite close 
to a central merger that happens in the plane of the sky, even 
though a small deviation from this geometry appears likely. 
For the sake of simplicity, we treat all other merger simula- 
tions in this paper as exactly central, with the merger axis 
aligned with the plane of the sky. 
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FIG. 9. — Comparison of the temperature profile across the shock front with 
deep Chandra observations of 1E0657-56. The sho ck-front is marke d with a 
vertical line. The crosses are the measurements by Markevitch 1 2006), with 
the leftmost point giving the temperature of the bullet. 



5.2. Baryon fractions of the clusters 

In our default models, we consider equal baryon fractions 
for the two merging clusters, given by ft, = 0.17, as sug- 
gested by r ecent cosmologica l parameter constraints based on 
WMAP3 dSpergel et al.ll2006l) . and we neglect the small frac- 
tion of baryons locked up in stars, so that all of the baryons 
are contained in the diffuse gas of the ICM. Clusters are gen- 
erally believed to be essentially closed-box systems due to 
their large mass, so the assumption of a baryon fraction equal 
to the universal value should be comparatively safe. On the 
other hand, there is some tension with respect to this in the 
low-density, low-era cosmological model favored by WMAP3, 
since detailed accounting of all observed baryonic compo- 
nents in_cJu£ters_jBndsJx>lead to values for f, smaller than 
0.17 dMcCarthv et al.ll2006h . 

However, in a number of test simulations where we varied 
f, =0.17, we found that our results are relatively insensitive 
to changing this parameter, unless we assign different baryon 
fraction to the clusters. The latter might be appropriate if the 
smaller sub-cluster has lost some of its baryons as a result of 
strong outflows in progenitor systems, or has locked up more 
of its baryons into collisionless stars. In this case we expect 
trends that are closely related to those found for changes in 
the concentration, which we considered earlier. If the baryon 
fraction of the bullet is lowered relative to that of the main 
cluster, it will be harder for the bullet's gas to pass through 
the main cluster's core. This should enlarge the separation 
between mass centroid of the bullet and the X-ray peak after 
core passage. Conversely, if the bullet is made more gas-rich 
while the parent's gas density is kept fixed or lowered, the de- 
celeration of the bullet by ram-pressure will become smaller, 
similar to the effect obtained when the concentration of the 
parent cluster is lowered, which will also tend to keep the 
mass centroid of the bullet and its X-ray cold front closer to- 
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FIG. 10. — The cumulative mass profile of our simulation model compared 
to gravitational lensing r esults for 1E0657- 56. The filled circles give the 
mass profile obtained by Bradac et al. 1 2006) for the parent cluster based on 
a combined weak and strong lensing analysis. The stars mark their inferred 
masses for the subcluster component. 



gether. 

We have carried out a number of additional simulations 
where we varied the baryon fractions in the two merging sys- 
tems, but otherwise adopted the parameters of our default 
model. The results of these simulations qualitatively con- 
firm the expectations described above without improving the 
match to the 1E0657-56 system in any significant way. For 
the sake of brevity, we therefore omit showing X-ray maps for 
these runs in a separate figure. 

5.3. Numerical resolution 

For all our default simulations we used a mass resolution of 
m dm = 6.225 x 10 8 M© and m gas = 1.275 x 10 8 M Q , translating 
to 4 x 10 6 particles for the parent cluster and 4 x 10 5 for the 
infalling bullet cluster. To test for any possible dependence on 
numerical resolution, we have rerun our default model with an 
increased particle number by a factor of 8, and an improved 
spatial resolution by a factor of 2 per dimension. This sim- 
ulation produces results for the strength of the shock and for 
the locations and velocities of the various hydrodynamic and 
dark matter features that are in excellent agreement with the 
simulation computed at our standard resolution. 

We conclude form this resolution test that the numerical 
resolution used in all our simulations is sufficient to provide 
numerically converged results for all the quantities studied 
here. 

5.4. Other parameters 

Other physical parameters that in principle should be ad- 
dressed in the simulations if one wants to find a 'per- 
fect' model for the 1E0657-56 system include the mass ra- 
tio of bullet and parent cluster, their intrinsic shape, the 
relative central entropy levels of the gas in the clusters, 
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FIG. 1 1. — The time evolution of total X-ray emissivity (top panel) and 
mean gas temperature (bottom panel) of our best-matching simulation model 
for 1E0657-56. The time corresponding to the observed moment is marked 
with a vertical dashed line, while the dotted lines give the core passages of 
the bullet. 



their rel ative baryon fractions, the effective viscosity of 
the gas flSji acki & Springel 2006), its thermal conductivity 
(Jubel gas et al.ll2004l) . a nd also non-thermal pressure compo- 
nents like cosmic rays ([P frommer et al.l 120061; lEnsslin et al.l 
2006; Jubelgas et al.l 120061) or intracluster magnetic fields 
(Dola g et al.ll2002l) . It would also be interesting to try to iden- 
tify a 1E0657-56 look-alike in a large-scale cosmological hy- 
drodynamical simulation, where the full dynamics of the in- 
fall and of gas accretion would be treated self-consistently. 

Obviously, this long list of possibilities quickly leads to a 
very large parameter space which will be difficult to constrain 
without additional observational constraints. Given that al- 
ready our simple toy model provides a good match to the 
key features that are observed in 1E0657-56, we refrain from 
studying the extended physical parameter space further in this 
paper, apart from one exception that we shall consider in the 
next section. 

6. MERGERS EST FIFTH FORCE MODELS 

Many extensions of the standard model of particle physics 
not only predict dark matter candidates in the form of new 
elementary particles, but also additional (long-range) fields 
which may couple diff erently to the dark matter particles and 
ordin a ry baryons (e.g. iDamour et aLlll990t iFarrar & Peebles! 
l2004t iGubser & Peeblesl 120041) . Interest in such models is 
strong also on astrophysical grounds. For example, they can 
produce cosmic voids that are 'empti er' of galaxies than stan- 
dard cold dark matter models (e.g. IFarrar & Peeblesl l2004t 




FIG. 12. — X-ray surface brightness profiles for a variation of our de- 
fault merger model where the encounter occurs with impact parameter b = 
12.5 kpc. The panel on top shows the system with a face-on orientation of 
the orbital plane. In the middle panel, the orbital plane is at 45 degrees with 
respect to the projection plane, in the bottom panel it is seen edge-on. 

iNusser et al.ll2005l) . which might be required to match obser- 
vational data. 

One of the simplest classes of such models can be phe- 
nomenologically characterized by a 'fifth forth' in the dark 
sector, mediated by a Yukawa-like potential 

^(r)=-/?^exp(-£) O) 

between dark matter particles of mass m, in addition to ordi- 
nary gravity. Here r s is an effective screening length of the 
scalar force, while (3 gives its strength on small scales rela- 
tive to ordinary gravity. The difficulty to explain a velocity 
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of t he bullet as high as ~ 4 700km s" 1 in ACDM cosmologies 
led lFarrar & Rosea (120061) to suggest that such a high speed, 
if true, may provide direct evidence for the existence of a fifth 
force. 

In this study, we have shown that the common approach 
of equating the shock velocity with the bullet velocity is in- 
correct, and that ACDM models in fact quite naturally pro- 
duce shock speeds in excess of 4500 kms -1 , despite the fact 
that the bullet is moving much slower than that, by about 
~ 2000km s" 1 . While this essentially eliminates the utility 
of 1E0657-56 to obtain strong constraints on the strength of a 
possible fifth forc e, we are here in t ereste d in explicitly testing 
the conjecture of Farrar & Rosenl (2006), i.e. we want to see 
by how much a fifth force can accelerate the bullet in direct 
simulations of the merging process. 

For definiteness, we choose a range of r s = 5.5 Mpc and ex- 
plore two different values of (3, namely (3 = 1.0 and 0.3. The 
size of r s is large enough to encompass both the parent clus- 
ter's virial radius and that of the infalling bullet cluster, such 
that the gravitational constant between dark matter particles is 
effectively boosted to (1 + j3)G when they are isolated. If we 
hence increase the random dark matter velocities in the dark 
halos of the initial clusters by a factor [((3+ 1 +/*,)/( 1 +fb)] 1 ^ 2 , 
we retain dynamical stability of the isolated systems in the 
fifth force model, with the same density profile as in our de- 
fault simulations. 

The range r s is also short enough such that a parabolic orbit 
from infinity leads only to a slightly increased velocity com- 
pared to a Newtonian infall when the virial radii first touch. 
We approximate this initial velocity of the encounter with the 
relative velocity two point masses would reach when falling 
together from infinity up to the initial relative separation of 
the sum of the two virial radii, which is the starting separa- 
tion adopted in our simulations. In doing so we implicitly 
assume that the relative binding forces between gas and dark 
matter hold the isolated systems together during the initial in- 
fall without appreciable shape distortions due to the violation 
of the equivalence principle in the fifth force model. We then 
expect a 17.2% increase in the relative velocity at the initial 
time for the (3 = 1 model, which we correspondingly include 
in our initial conditions for the fifth force simulation. We note 
that for an infinite range of the fifth force, this correction fac- 
tor would increase to 30%. 

In FigureQj] we show maps of the X-ray surface brightness 
with overlaid mass contours for the simulations with j3= 1.0 
and (3 = 0.3, respectively, at a time when the mass centroids 
of the bullet and the parent cluster are 720 kpc apart. The X- 
ray morphology seen in the maps can be directly compared to 
the corresponding maps for ordinary mergers shown in Fig. [7] 
What is immediately apparent is the rather large offset be- 
tween the dark matter and the gas of the bullet in the j3 = 1.0 
case. The centroid of the dark mass of the bullet is 230kpc 
ahead of the shock front in this case, and 300 kpc ahead of the 
edge of the bullet. An analysis of the bow shock shows that 
the Mach number of the shock is slightly lower in this fifth 
force run, with M. = 2.5, but the preshock temperature is a 
bit higher (10 keV), such that the inferred shock velocity is 
~ 3900km s" 1 . Also in this case there is a preshock inflow ve- 
locity, but with a slightly reduced size of ~ 700 km s" 1 , which 
in part is explained by the earlier time of the 'best-match', 
here defined in terms of the separation of the mass centroids 
of bullet and parent cluster. 

The different timing of the merger is shown in more detail 
in Figure [141 where we plot the coordinates of the mass cen- 




FlG. 13. — Map of the X-ray emission of a merger model with a 5th-force 
in the dark matter sector, characterized by j3 = 1.0 (top panel) and /3 = 0.3 
(bottom panel). The projected mass distribution is given by overlaid contours. 
In both cases, the simulations are shown at a time when the two mass peaks 
are separated by 720 kpc. 



troids of bullet and parent cluster as a function of time, com- 
paring results for j3 = 1, (3 = 0.3 and (3 = 0) (i.e. no fifth force). 
The times when the separations reach 720kpc are shown with 
vertical dotted lines. Clearly, the (3 = 1.0 model leads to a 
significantly earlier core passage and matching time than the 
other models, as a result of the higher velocity reached by 
the bullet. In fact, at the time of best-match, the dark mat- 
ter of the bullet is moving with ~ 3800 kms -1 , which is here 
by chance quite close to the numerical value obtained for the 
shock speed. This corresponds to a substantial increase of 
Avboost = 1200 kms" 1 relative to the value predicted for the 
corresponding ACDM run. If (3 is lowered to (3 = 0.3, the 
velocity increase is reduced to Avboost = 410kms -1 , and the 
opening angle of the bow shock becomes larger (see Fig.[T3l. 
The latter is hence in principle an interesting diagnostic for 
the strength of the fifth force. 

The above clearly shows the substantial impact models with 
non-trivial interactions in the dark sector have on cosmolog- 
ical structure formation, and in particular on the non-linear 
regime that governs the internal structure and dynamics of ha- 
los. So far, comparatively little simulation work has addressed 
this problem (with notable exceptions, e.g. lNusser et al.ll2005l: 
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FIG. 14. — Positions of the dark matter mass peaks of bullet and parent 
cluster as a function of time. We compare a 5th-force model with (3 = 1.0 
(solid lines) and (3 = 0.3 (dot-dashed lines) with our default model (dashed 
lines). The times where the separation of bullet and parent cluster are equal to 
720 kpc are marked with dotted vertical lines. For the fifth force models, the 
bullet reaches a higher velocity, and core passage happens earlier. When the 
separation is 720 kpc, the bullet's dark matter moves with 3800 kms~' in the 
j3 = 1.0 case and with 3010km s for /3 = 0.3, while without the additional 
scalar force, the velocity reaches only 2600kms~' at a time with an equal 
separation between bullet and parent cluster. 



iKesden & Kamionko wski 2006), but as seen here, direct nu- 
merical simulation methods are a powerful tool to study the 
implications of such models and to constrain them with ob- 
servational data, so future work in this area should be very 
worthwhile. 

7. CONCLUSIONS 

The bullet cluster 1E0657-56 offers a unique perspective 
into cluster physics and the hierarchical build-up of cosmic 
structure. Deep X-ray observations have revealed a strong 
bow shock in front of the bullet, and a contact discontinu- 
ity, a 'cold front', that gives its bright central region a sharp, 
wedge-like edge. In addition, mass reconstructions based on 
weak- and strong gravitational lensing have allowed a reli- 
able determination of the projected mass distribution of the 
system. The latter features a spatial offset between the mass 
centroids and the baryonic mass centroids, which represents 
strong evidence for the existence of a non-baryonic dark mat- 
ter component. 

This combination of powerful observations destines 
1E0657-56 to serve as laboratory for constraining funda- 
mental aspects of cluster physics and of our cosmological 
model. For example, the system has already been used to 
determ ine constraints on the co llisional cross-section of dark 
matter (Markev itch et al.ll2004l) . It also has become one of the 
hardest challenges for alternative gravity models like MOND 
that try to circumvent the need for dark matter. While claims 
have been made that MONDian-like th eories could possi- 
bly reproduce the bullet cl uster as well (|Angus et alj 120061; 
Brownstein & Moffat 2007), nothing that comes close to the 



quality of the matches presented here for ACDM has been 
presented thus far to support such scenarios. 

The point we have focused on in this paper is the high in- 
ferred shock velocity in the bullet cluster, which is a direct 
consequence of the measured Mach number and the pre-shock 
temperature of the gas. The literature thus far has assumed 
that this shock velocity can be directly identified with the ve- 
locity of the bullet subcluster itself. Theoretical papers have 
then pointed out that such a large velocity for the bullet is un- 
comfortably Jri^h 2 _ejtherjm2lyirm the system is a rather 
rare outlier ( Ha vashi & Whitd2006|) or may need new p hysics 
in the dark sector to be explained dFarrar & Rosen 200_6|). 

Using explicit hydrodynamical simulation models of clus- 
ter mergers, in 3D and with self-consistent gravitational dy- 
namics of dark matter and baryons, we have here shown that 
a direct correspondence between the inferred shock velocity 
of the bullet and that of its associated mass is however not 
given. In the case of 1E0657-56, the two velocities can dif- 
fer both because the upstream gas is accelerated towards the 
shock-front by gravitational effects that can be felt ahead of 
the hydrodynamical shock, and secondly, because the shock 
front is moving with a different speed than the mass cen- 
troid of the bullet subcluster. Our simulation model based 
on a 1:10 merger of ordinary ACDM halos on a central, zero- 
energy orbit reproduces the inferred shock speed well, even 
though the subcluster moves only with 2600km s" 1 . At the 
same time, the model provides surprisingly good matches to 
the observed temperature profile across the shock and cold 
fronts, and to the observational determined mass distributions 
in the gas from X-rays, and in the total mass from lensing. We 
conclude that the bullet cluster 1E0657-56 is consistent with 
the simplest possible scenario one can construct in ACDM for 
a 1 : 10 mass-ratio merger between a very rich cluster of galax- 
ies and a smaller sub-clus ter. In fact, based on the results 
of lHavashi & White! (120061) we estimate that the likelihood to 
find a most massive dark matter substructure with a veloc- 
ity of at least 2860 km s -1 embedded in a parent system with 
y 200 = 1770kms _1 at z = 0.296 is as high as / ~ 7%. Here we 
conservatively accounted for the increase of the relative ve- 
locity of the subcluster relative to the parent cluster, which is 
bounded to be less than 260 kms -1 , depending on how much 
mass the subcluster has already lost. 

We have also shown that the size of the spatial offset be- 
tween mass and baryonic peak is quite sensitive to the struc- 
tural details of the merging systems. While a spatial offset 
between the mass centroids of dark matter and baryons arises 
naturally as a result of the hydrodynamic ram pressure ex- 
perienced by the bullet's gas during core passage, this dis- 
placement depends on the relative gas fractions and concen- 
trations of the merging clusters, in addition to being highly 
time-variable. This could also be the reason why some ear- 
lier numerical model s of the bullet c luster missed to detect 
the offset altogether dTakizawal 120051) . As Figure [4] shows, 
for low initial concentrations of the parent cluster, the shock 
front and the edge of the bullet can even overtake the mass 
centroid of the cluster at late times during the mer ger, thereby 
revers ing the sign of the displacement seen by IClowe et all 
(2006). One is hence not necessarily guaranteed to see a spa- 
tial offset of similar size as in 1E0657-56 between dark mass 
centroid and gaseous centroid in all massive cluster mergers, 
and during the whole post core-passage phase. Our results 
emphasize that an interpretation of highly non-linear systems 
such as 1E0657-56 requires self-consistent simulation mod- 
els that account accurately both for the gravitational and hy- 
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drodynamical dynamics of dark matter and baryons during 
cluster mergers. 
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